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ABSTRACT 


MODIS collects data in both the reflected solar and thermal emissive regions using 36 spectral bands. The center 
wavelengths of these bands cover the3.7 to 14.24 micron region. In addition to using its on-board calibrators (OBC), 
which include a full aperture solar diffuser (SD) and a blackbody (BB), lunar observations have been scheduled on a 
regular basis to support both Terra and Aqua MODIS on-orbit calibration and characterization. This paper provides an 
overview of MODIS lunar observations and their applications for the reflective solar bands (RSB) and thermal emissive 
bands (TEB) with an emphasis on potential calibration improvements of MODIS band 21 at 3.96 microns. This spectral 
band has detectors set with low gains to enable fire detection. Methodologies are proposed and examined on the use of 
lunar observations for the band 21 calibration. Also presented in this paper are preliminary results derived from Terra 
MODIS lunar observations and remaining challenging issues. 
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1. INTRODUCTION 

Since launch in December 1999, more than 10 years of global earth observations have been made by Terra MODIS and 
used to derive a wide range of data products for studies of changes in the Earth’s system of land, oceans, and 
atmosphere. MODIS is a multispectral filter radiometer and collects data in 36 spectral bands. Twenty of the spectral 
bands are in the reflected solar region and the other 16 are in the thermal emissive region. The reflective solar bands 
(RSB), covering wavelengths from 0.41 to 2.2pm, are calibrated on-orbit using a solar diffuser (SD) and solar diffuser 
stability monitor (SDSM). The thermal emissive bands (TEB), with wavelengths spanning 3.7 to 14.4pm, are calibrated 
by a blackbody (BB). In addition to its on-board calibrators (OBC), MODIS lunar observations, which are scheduled on 
a near monthly basis, are used to support its on-orbit calibration and characterization 1 " 3 . Figure 1 shows the MODIS scan 
cavity and OBC. MODIS also has a space view (SV) port designed for observations to correct instrument background 
and offset. Regularly scheduled lunar observations are also made through the SV port. 

MODIS TEB detectors are calibrated using their responses to the on-board BB on a scan-by-scan (1.4771 sec) basis, 
with an exception of band 21 at 3.96pm. This band was designed for fire detection with a specified maximum 
temperature of 500K. At this temperature, the band 21 spectral radiance is more than 2 orders of magnitude higher than 
that from the BB nominal operating temperature of 290K for Terra MODIS. The small detector responses, or digital 
counts, of band 2 1 to the on-orbit calibration source (BB) coupled with their low signal-to-noise ratios (SNR) are 
inadequate to derive accurate and stable calibration coefficients for band 21 over its entire dynamic range. 


Consequently, band 2 1 detector calibration coefficients (or gains) are derived using data from periodic BB warm-up and 
cool-down operations during which the BB temperatures vary from instrument ambient of approximately 272K to 3 15K. 
These calibration coefficients are detector and mirror side dependent 4 " 5 . Figure 2 illustrates Terra MODIS bands 21 
(3.96pm), 22 (3.96pm), 31 (11pm), and 32 (12pm) middle detector responses in digital numbers (dn BB ) as a function of 
calibration radiances (AL BB ) with data from the BB cool-down operation on day 2004115. The calibration radiances 
AL bb have included contributions from the BB, scan mirror, and instrument cavity. The corresponding detector 
responses dn BB are corrected for instrument background provided by the SV observations. Even at 315K, band 21 
detector responses are still very small (less than 100 counts) compared to its dynamic range upper limit. Because of this, 
MODIS band 2 1 is calibrated using a simple linear algorithm. 



Figure 1: MODIS internal scan cavity, its on-board calibrators (OBC), and lunar observations through the space 
view (SV) port. 
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Figure 2: Terra MODIS bands 21, 22, 31, and 32 detector responses in digital numbers (dn BB ) as a function of 
calibration radiances (AL BB ). Data from BB cool-down on day 2004115. 

Except for band 21, the dynamic ranges of bands 22, 31, and 32 are suitable for deriving high quality calibration 
coefficients. Based on the design requirements, band 22 has the same relative spectral response (RSR) as band 21 but 



with much higher detector gains. As shown in Figure 2, band 22 detector gain is approximately 25 times higher than 
band 21. 

This paper provides an overview of MODIS lunar observations and their applications in support of both the RSB and 
TEB on-orbit calibration and characterization, including radiometric calibration stability monitoring, spatial 
characterization, and crosstalk (optical and electronic) assessments. In particular, this paper discusses potential 
calibration improvements for MODIS band 21 through the use of its lunar observations. Preliminary results derived 
from Terra MODIS lunar observations are presented. Remaining challenging issues and potential improvements for 
future studies are also discussed. The approaches developed from MODIS lunar calibration can also be applied to other 
sensors, such as the NPOESS Preparatory Project (NPP) and Joint Polar Satellite System (JPSS) VIIRS. 


2. OVERVIEW OF MODIS LUNAR OBSERVATIONS 

Since lunar surface reflectance properties are known to be extremely stable, MODIS lunar observations were planned 
initially to monitor the RSB on-orbit radiometric calibration stability 3 . A special lunar planner tool was developed by the 
MODIS Characterization Support Team (MCST) in support of its lunar calibration activities. Since launch, lunar 
observations have been scheduled and implemented on a regular basis for both Terra and Aqua MODIS. Because the 
lunar surface irradiance strongly depends on the viewing geometry, the scheduled lunar observations have been made at 
the same phase angle, 55° to 56° (waning) for Terra MODIS and -55° to -56° (waxing) for Aqua MODIS. In order to 
meet the phase angle requirement, spacecraft maneuvers (rolls) are often needed. MODIS lunar observations are made 
through its SV port with data sector rotation operations performed right before and after the lunar observations. 

Figure 3 shows examples of Terra MODIS band 1 (0.65 pm) lunar observations made during a period from September to 
December 2006. Thus far, more than 95 and 70 lunar calibration events have been scheduled and carried out for Terra 
and Aqua MODIS, respectively. In addition, both MODIS instruments have acquired many un-scheduled lunar 
observations. In general, the unscheduled lunar observations without spacecraft maneuvers have different phase angles. 
After corrections for the view geometry differences, the unscheduled lunar observations could also be used for the RSB 
radiometric stability monitoring and other applications 6 . 
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Figure 3: Terra MODIS scheduled lunar observations (band 1 at 0.65pm) from September to December 2006 


3. APPLICATIONS OF LUNAR OBSERVATIONS 


3. 1 Lunar Observations for RSB 

MODIS RSB on-orbit calibration is performed using its on-board SD/SDSM system. The lunar observations are also 
made to monitor the RSB radiometric calibration stability. MODIS views the sunlight reflected from the SD and the 
lunar surface via the same optical system. The only difference is that the SD and lunar observations are made at two 
different angles of incidence to the scan mirror. Because of this, MODIS lunar observations are used together with SD 
observations to determine changes in RSB response versus scan angle (RVS) 7 " 8 . The SD calibration coefficients, (m^so, 






are derived from detector responses to the SD with corrections for instrument background, temperature, and a number of 
view geometry-dependent factors. Similarly, lunar calibration coefficients, (mi) Moon , require corrections for the lunar 
view geometry. The RSB calibration coefficients are derived for each band, detector, sub-sample, and mirror side using 
the following expressions, 


(mi)sd ~ F S d/ (dn*) S D, ( 1 ) 

(wi l) Moon F][f oon / (dfl )mooh • ( 2 ) 

where F SD is a correction factor, which depends on the SD BRF, thus the SD viewing geometry, its on-orbit degradation, 
the SD screen transmission if it is in place, and the Earth-Sun distance. F Moon , is a correction factor that depends on the 
lunar surface reflectance properties, which are viewing geometry dependent, the Earth-Moon distance, and the Moon- 
Sun distance. (dn*) SD and (dn*) Moon are the detector digital responses to the SD and Moon, respectively, with instrument 
background, viewing angle, and temperature effects corrected. To minimize the effect due to different lunar viewing 
geometry, lunar observations have been carefully scheduled with the same phase angle for each MODIS instrument. 

MODIS has 36 spectral bands and a total of 490 individual detectors. They are located on four focal plane assemblies 
(FPA) as shown in Figure 4. All MODIS spectral bands are aligned in the cross-track direction and detectors of each 
band are aligned in the along-track direction. When MODIS views the Moon, the same lunar image gradually moves 
across all the bands (detectors) on the FPA. The normalized lunar profiles from individual detectors in both directions 
can be used to derive the detector-to-detector registrations (DDR) and band-to-band registration (BBR). Recently, 
MOST has developed an approach to use lunar observations to characterize the modulation transfer function (MTF) for 
the RSB. In addition, MODIS lunar observations have been used for on-orbit characterization of its SWIR detector 
thermal leak and electronic crosstalk 9 . 



Another useful application of lunar observations is to enable sensor calibration inter-comparison 10 . This approach can 
examine calibration differences among sensors and, therefore, help understand and resolve calibration discrepancies 
between sensors. It can consistently track sensors’ radiometric stability even if they are operated on different platforms 
or over different timeframe. For radiometric stability monitoring and calibration inter-comparison purposes, a lunar 
radiometric model, such as the one from USGS is needed to provide accurate and consistent corrections for the lunar 
viewing geometry 11 . Results from MODIS lunar observations show that Terra and Aqua MODIS have been consistently 




calibrated to within ±1% for most RSB. The BBR results derived from MODIS lunar observations generally agree with 
that derived from its on-board device, a spectroradiometric calibration assembly (SRC A). 

3.2 Lunar Observations for TEB 

With the same philosophy, lunar observations can be applied for the TEB calibration stability monitoring. However, 
earth-observing sensors normally saturate when viewing the Moon since typical earth scene brightness temperatures are 
much lower than the Moon. Table 1 is a summary of MODID TEB key design parameters, including their center 
wavelengths and bandwidths, typical and maximum specified radiances and temperatures, and noise equivalent radiance 
and temperature differences. Only band 21, designed for fire detection, has specified maximum temperature above 
335K. For Terra MODIS, however, the actual saturation temperatures, also listed in Tablet, are approximately 392, 387, 
and 374K for bands 31, 32, and 36, respectively. 


Table 1: Terra MODIS TEB specifications (CW: center wavelength, BW: bandwidths, Ltyp/Lmax: typical/maximum 
specified radiances, T Ltyp/T Lmax: temperature at typical/maximum radiance, NEdL/NEdT: noise equivalent 
radiance/temperature difference at typical radiance, T Sat: sensor actual saturation temperature) 


Band 

CW 

(ran) 

BW 

(ran) 

Ltyp Lmax 

(W/m A 2/sr/D) (W/m A 2/sr/D) 

NEdT 

(K) 

NEdL 

(W/m A 2/sr/D) 

TLtyp 

(K) 

T Lmax 
(K) 

T Sat 
(K) 

20 

3750 

180 

0.45 

1.71 

0.05 

0.0010 

300 

335 

335 

21 

3959 

59 

2.38 

85.44 

0.20 

0.0154 

335 

500 

478 

22 

3959 

59 

0.67 

1.89 

0.07 

0.0019 

300 

328 

329 

23 

4050 

61 

0.79 

2.16 

0.07 

0.0022 

300 

328 

330 

24 

4465 

65 

0.17 

0.34 

0.25 

0.0022 

250 

264 

317 

25 

4515 

67 

0.59 

0.88 

0.25 

0.0062 

275 

285 

316 

27 

6715 

360 

1.16 

3.21 

0.25 

0.0108 

240 

271 

323 

28 

7325 

300 

2.19 

4.47 

0.25 

0.0172 

250 

275 

319 

29 

8550 

300 

9.59 

14.55 

0.05 

0.0090 

300 

324 

330 

30 

9730 

300 

3.70 

6.34 

0.25 

0.0219 

250 

275 

358 

31 

11030 

500 

9.56 

13.26 

0.05 

0.0070 

300 

324 

392 

32 

12020 

500 

8.95 

12.10 

0.05 

0.0061 

300 

324 

387 

33 

13335 

300 

4.53 

6.56 

0.25 

0.0183 

260 

285 

334 

34 

13635 

300 

3.77 

5.03 

0.25 

0.0161 

250 

268 

341 

35 

13935 

300 

3.11 

4.42 

0.25 

0.0141 

240 

261 

341 

36 

14235 

300 

2.08 

2.96 

0.35 

0.0154 

220 

238 

374 


Since not all the lunar pixels have the same temperatures and the retrieved brightness temperatures vary significantly 
across the entire lunar surface, only a fixed number of pixels with the highest temperatures are used to track each 
detector’s long-term change caused by its on-orbit calibration. This approach has been applied to Terra MODIS bands 
31 and 32. Results from 10 years of lunar observations indicate that Terra MODIS TEB calibration has been generally 
stable with no obvious drift over its entire mission. For Aqua MODIS, the saturation temperatures are much lower than 
Terra MODIS and, therefore, this approach cannot be directly applied to Aqua MODIS. 

For TEB bands with saturated lunar pixels, we plan to use a fixed number of pixels in the middle temperature range to 
monitor the changes in detector responses over time. A ratioing approach has been applied to some reflective solar 
bands that have saturated lunar pixels. A similar approach can also be applied to the TEB in order to characterize their 
detector-to-detector calibration differences. Since most TEB bands in both Terra and Aqua MODIS saturate when 
viewing the Moon, the lunar BBR characterization is limited to their VIS and NIR spectral bands. 

It was identified from pre-launch that Terra MODIS had experienced an optical leak among TEB photoconductive (PC) 
detectors (bands 31-36). On-orbit lunar observations have been routinely used to derive the optical leak correction 
coefficients applied to the LIB algorithm 9 . 



4. MODIS BAND 21 CALIBRATION 


4. 1 Methodology 

Except for band 21 which applies a simple linear algorithm, MODIS TEB calibration uses a quadratic algorithm with 
reference to the on-board BB on a scan-by-scan basis. As illustrated in Figure 2, band 21 detector responses to the BB 
and associated SNR are too low to produce useful calibration parameters that can be used to accurately retrieve the high 
temperature scenes. In this study, an approach is developed as an attempt to apply lunar observations to improve B21 
calibration over its entire dynamic range. 

For a quadratic algorithm, calibration source radiance (space background and instrument thermal emission subtracted), 
AL S , is related to the corresponding detector responses (background responses subtracted), dn s , through a set of 
quadratic calibration coefficients (a 0 , bi, and, a 2 ), 

AL S = a 0 + b,dn s + a 2 (dn s ) 2 (3) 

Equation 3 can be applied to on-board BB and lunar observation. For TEB on-orbit calibration, AL BB includes three 
terms with contributions from the on-board BB, scan mirror, and instrument cavity. For lunar observations, AL Moon has 
only two components: one is due to the thermal emission from the Moon and another due to solar reflection from the 
lunar surface. Thus 

ALMoon = e M oonRVSsvL(T M oon, + solar_reflection (A) (4) 

where 

ZMoon is the lunar surface emissivity at wavelength A, 

RVS S v is the response versus scan angle (RVS) at the sensor’s space view (SV) angle, 

Tmoou is the temperature of the Moon on a pixel-by-pixel level 

In general, the reflected solar contribution is extremely small for the longer wavelength TEB bands, such as band 3 1 . It 
needs, however, to be considered for the bands at relatively short wavelengths, such as band 21. In order to use 
Equations 3 and 4 for band 21 calibration, the following steps should be implemented: 

1 . Apply Equations 3 and 4 to band 3 1 and compute T Moon for each pixel 

2. Derive lunar emissivity for band 21 wavelength 

3. Derive solar reflection for band 21 wavelength 

4. Apply Equations 3 and 4 for band 21 calibration 

The first step is straight forward as the lunar emissivity at band 31 wavelength can be found from the literature and its 
solar reflectance term can be ignored 12 " 14 . To derive the lunar emissivity at the wavelength of band 21, a best fit between 
band 31 and band 21 spectral radiances from the moon (which depend on lunar emissivity) and to that from the on- 
board BB is performed. After steps 1 and 2, the solar reflection at band 21 wavelength can be determined at the 
relatively high end of its detector responses. 

4.2 Preliminary Results 

Figure 5 shows several curves of Terra MODIS band 21 radiances versus band 31 radiances (the right side plot is a 
zoom-in of the left side) for their detector 5 (production order). In this fitting, 0.9 is used as the lunar emissivity for band 
3 1 and its solar reflection is considered to be negligible. The thick red curve represents data from BB warm-up and cool- 
down in 2002. The other four curves are constructed with different lunar emissivity and solar reflection values for band 
21. In Figure 5, the best lunar emissivity is estimated to be about 0.682. Although the solar reflection is included as a 
fitting variable, it is still found to be zero after obtaining the best matching of band 2 1 and band 3 1 lunar radiances with 
their BB radiances. 
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Figure 5: Terra MODIS band 21 and band 3 1 (mirror side 5 and detector 5) BB and lunar radiances (see text for 
details). Radiance unit: W/m 2 /sr/p. 


An alternative approach is to combine steps 2 and 3 and determine the lunar emissivity and solar reflection at the same 
time by best fitting band 21 lunar data to a curve (line) projected based on band 21 BB calibration data. Two examples 
are presented in Figure 6 with different lower radiance limits (LR). In this analysis, the emissivity of band 31 is again 
set to be 0.9. BB WUCD data in 2008 is used. Only detector 5 results are illustrated here. The best emissivity and 
reflectance estimates are [0.686, 1.169] and [0.674, 1.261] with LR=2 and LR=4, respectively. This study has examined 
different cases by using different emissivity for band 3 1 and different lower radiance limits for band 2 1 . Results from 8 
out of 10 band 21 detectors (two end detectors excluded) are provided in Figure 7 to help assess the methodology and 
overall impact due to different input. 

Terra (msl) D-5 Terra (msl) D-5 
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Figure 6: Band 21 lunar emissivity and solar reflection derived with different radiance ranges. Left: lunar data 
with lower radiances above 2 (W/m 2 /sr/p) are used for best fitting; Right: lunar data with lower radiances above 4 
(W/m 2 /sr/ju) are used for best fitting. 


Some small but noticeable differences among individual detectors are noticed in their retrieved parameters. This is 
likely due to small differences of individual detectors’ view geometry, small time and over-sampling factor differences, 
and, of course, their calibration differences or uncertainties. With the increase of lower radiance limits, the detector-to- 
detector differences become much larger. Considering all the features in Figure 7, one could select the best band 21 
parameters required for implementing its lunar calibration. With different band 31 emissivity values, the changes in the 


band 21 solar reflection term are generally small. On the other hand, the best estimated emissivity is strongly dependent 
on the knowledge (quality) of band 31 emissivity. 




Radiance Low-Level (B21) Radiance Low-Level (B21) 

Figure 7: Terra MODIS band 21 (8 middle detectors) lunar emissivity (left) and solar reflection (right) estimations with 
different band 3 1 emissivity values and different best fitting lower radiance limits. Reflectance (label) above refers to 
the solar reflection (Eqn. 4) contribution to the total of lunar radiance of band 21 (in W/m 2 /sr/p). 


To improve the band 21 emissivity estimation, we have applied the same approach used in Figure 5 to band 22. This is 
based on the consideration that band 22 BB calibrations are of higher quality and more stable than band 2 1 while both 
bands have nearly the same relative spectral response (RSR). The results from this approach are presented in Figure 8. 
Again the solar reflection is found to be zero after performing the best matching of band 22 and band 31 lunar radiances 
to their BB radiances. The best lunar emissivity estimation for band 22 is about 0.681, very close to the band 21 value 
shown in Figure 5. Another improvement planned for future study is to separate the lunar data with reflected sunlight 
from those without. For the data without reflected sunlight, we can ignore the solar reflection contribution and focus on 
the emissivity estimation. 
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Figure 8: Terra MODIS band 22 and band 31 BB and lunar radiances (see text for details). Radiance unit: 
W/m 2 /sr/p. 


5. SUMMARY 


Lunar observations have been successfully used in support of MODIS on-orbit calibration and characterization for both 
the reflective solar bands (RSB) and thermal emissive bands (TEB). Applications include using the Moon for sensor 
radiometric calibration stability monitoring, calibration inter-comparison, band-to-band registrations (BBR) and MTF 
characterization, and sensor crosstalk (optical and electronic) evaluation. This paper has provided an overview of 
MODIS lunar observations and their applications developed by the MODIS Characterization Support Team (MCST). It 
has discussed potential improvements of band 21 on-orbit calibration quality at high temperatures. Different algorithms 
to estimate band 21 lunar emissivity and solar reflection have been proposed and examined using Terra MODIS lunar 
observations. Preliminary results indicate that, with additional effort and data analysis improvements, the approaches 
developed in this paper can lead to good estimates of band 21’s lunar emissivity and solar reflection, and thus enable the 
use of lunar observation to improve band 2 1 calibration. 
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